arXiv:1502.07344vl [astro-ph.GA] 25 Feb 2015 


Mon. Not. R. Astron. Soc. 000 , [l][13] (2009) Printed 27 February 2015 (MN IAI^X style file v2.2) 


Galaxy formation with radiative and chemical feedback 


L. Graziani 1 *, S. Salvadori 2 , R. Schneider 1 , D. Kawata 3 , M. de Bennassuti 1 , A. Maselli 4 

1 INAF Osservatorio Astronomico di Roma, Via Frascati 33, 00040, Monte Porzio Catone (RM), Italy 

2 Kapteyn Astronomical Institute, Landleven 12, NL-9747 AD Groningen, the Netherlands 

3 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK 

4 EVENT Lab for Neuroscience and Technology, Universitat de Barcelona, Passeig de la Vail d’Hebron 171, 08035 Barcelona, Spain 


April 2014 


ABSTRACT 

Here we introduce GAMESH, a novel pipeline which implements self-consistent ra¬ 
diative and chemical feedback in a computational model of galaxy formation. By com¬ 
bining the cosmological chemical-evolution model GAMETE with the radiative transfer 
code CRASH, GAMESH can post process realistic outputs of a N-body simulation de¬ 
scribing the redshift evolution of the forming galaxy. After introducing the GAMESH 
implementation and its features, we apply the code to a low-resolution N-body sim¬ 
ulation of the Milky Way formation and we investigate the combined effects of self- 
consistent radiative and chemical feedback. Many physical properties, which can be 
directly compared with observations in the Galaxy and its surrounding satellites, are 
predicted by the code along the merger-tree assembly. The resulting redshift evolution 
of the Local Group star formation rates, reionisation and metal enrichment along with 
the predicted Metallicity Distribution Function of halo stars are critically compared 
with observations. We discuss the merits and limitations of the first release of GAMESH, 
also opening new directions to a full implementation of feedback processes in galaxy 
formation models by combining semi-analytic and numerical methods. 

Key words: Cosmology: theory, galaxies: formation, evolution, stellar content, star: 
formation, Population II, reionisation, radiative feedback, Milky Way 


1 INTRODUCTION 


Since the pioneering work of ICouchman fe Reel d 1986h . 
the radiative feedback induced by cosmic reionisation 
has been recognised to have a strong impa ct on 


the formation and ev olution of galax i es (IGn cdin 
Benson et al.l l2002bllal: [Somervilldl l2002i; iHoeft et al. 


J] 

Okamoto^^oj^TheurRJgOO^I Salyadori fe Ferrara! 2009 : 

Lunnan et al.ll2012l : ISawala et al.ll2014l '). imprinting specific 
signatures in the observed properties of the Local Universe 
llBrown et al.l[2oTd ). 

Semi-analytic mode ls of galaxy formation (SAM) (see 
e.g. iBenson et al . 1 l2002al lbl or ISomervilB l2002j l and numer- 


2000 : 


2006 ; 


ical s im ulations, both un i- dimen sio nal ( Thoul_J^WeinbeG 


19961: iKitavama et all 


Gnedinl 200d: Hoeft 


et al. 


200d : iDiikstra et al.l 120041 : 


Sobacchi fc Mesingeil 1201 .i) and three-dimensional (see 


2006i ; lOkamoto, Gao fe Theunsl 


20081 and references therein) have been performed to 


determine the efficiency of the radiative feedback and 
the way it affects galactic star formation (SF) by either 
an extended ( Gradual ) or instant (Sharp) process in red- 


* E-mail: luca.graziani@oa-roma.inaf.it 


shift; ho wever a physically mo tivated prescription is still 
missing (jNoh fe McQuinnll2014l 'l. Models with a sharp SF 
suppression after reionisation generally fail to reproduce 
the present, observed luminosity function and predict 
an unacceptable gap between the faint and the bright 
populations, also showing that a redshift dependence in 
modelling feedback effects is necessary. The value of the 
critical mass below which galaxies should be strongly 
affected by photo-ionisation at fixed redshift is also subject 
to an intense debate: recent numeric al studies (IHoeft et al.l 
I2006S : lOkamoto. Gao fc Theunsll200 8h find about one order 
of magnit ude lowe r masses compared to the original cal¬ 
culations llGnedinll2000h . but all the conclusions critically 
depend on the reliability of the ionising background (UVB) 
adopted by the SAM. 


The effects of an UVB maintaining the Universe re¬ 
ionised can be also considerable on the structure and prop¬ 
erties of the Local Group. The tendency of numerical sim¬ 
ulations to over-predict the substructures on small scales 


ing satellite problem" (Moore et al. 199S 

: Klypin et al. 1999:; 

Bovlan-Kolchin, Bullock & Kanlinghat 

2013 ), is often in- 


terpreted as inability to reproduce the inefficient or sup- 
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pressed star formation acting in the group of local dwarf 
galaxies, causing small h alos to remain "dark", and then 
invisible to observations (lAlvarez et al.l 120091 : iBusha et al.1 
l20ld : lOcvirk et al.ll2014| b Feedback processes have been in¬ 
voked to j ustify this f ailu re of the ACDM model, both me- 


2009[ 201ol) and radiative (Efstathiou 1992: Gnedin 

2000; 

Benson et al. 2002al; Somerville 20Cbf Madau et al.l 

2008; 

Noh & McQuinm 2014: Milosavlievic & Bromm 20141). 


A correct determination of the epoch of hydrogen reioni¬ 
sation on the scale of galaxy formation, self-consistently with 
the SF suppression operated by radiative feedback, could 
shed some light also on the nature of the ultra-faint dwarfs 
population found in t he Sloan Digital Sky Survey (SPSS, 


Willm an et al.1 120051 : iBelokurov et al.l l2006t IZucker el, al.l 
2006bl al. The discovery of these galaxies, roughly doubling 


the number of known Milky Way (MW) satellites, forced 
an update of the old models and opened new questions on 
whether this population is made by pre-reionisation fossils 
or it is just the lowest-luminosity tail of classical dwarfs. 
The answer will critically depend, again, on the way the gas 
ionisation is modelled and on how radiative effects are taken 
into account. 

By using a data-calibrated model for the forma¬ 
tion of the MW and its dwarf satellites, which includes 
the presence of inefficient star-forming mini-halos and a 
heuristic prescription to account for radiative feedback, 
ISalvadori fe Ferraral (l2009l l have been able to simultane¬ 
ously reproduce the observed iron-luminosity relation and 
the Metallicity Distribution Function (MDF) of nearby 
dwarf galaxies, including the ultra-faint population. The 
authors argued that ultra-faint dwarf galaxies might rep¬ 
resent the fossil relics of a once ubiquitous population 
of H 2 -cooling mini-halos that formed before reionisation 
(z > 8.5), simil arly to what has been found by inde¬ 
pendent groups jMadat^et^fJJ^OO^jKoi 3QSOV et al.| | 2009l: 


__ lOP 

Boyill fc Ricottil 120091 : Munoz et al.l 120091 Bovill fc Ricottil 
2011allbl) and in agreement with a series of more recent 


studies llSalvadori fc Ferraral 120121 ; ISalvadori et al.l I20l4) . 
In particular, they find that a gradual suppression of the 
star-formation in increasingly massive (mini-)halos is nec¬ 
essary to match the observed iron-luminosity relation and 
the MDF of nearby dwarfs. Although this empirical re¬ 
lation is consistent with the minimum mass for SF set- 
tled b y an a dvanced s emi-analytic treatment of reionisation 
llSalvadori et al.ll2014r ). the results still rely on the assump¬ 
tions made on the nature and efficiency of feedback acting 
on SF halos. A more detailed modelling is certainly required 
to understand the local effects of in-homogeneous radiative 
feedback, and its impact on the luminosity function of MW 
satellites. 

It should be noted, on the other hand, that the ob¬ 
served luminosity function of MW satellites is certainly not 
sufficient to fully constrain semi-analytic m odels as shown, 
for instance, by alternative approaches llMaccid et al.l [2OI0I : 
ILL De Lucia fc Helm ill 201 (J l. The observed number and dis¬ 
tribution of luminous satellites can be reproduced with ac¬ 
curacy also by assuming a sharp (i.e. instant) reionisation 
and a strongly mass-dependent star formation efficiency or 
by invoking other physical processes as supernovae feedback 
or tidal stripping._ 

Recent investigations llAlvarez et alj|2009t IBusha et al.l 


120101 : iLunnan et al.ll2012l : iLi et al.ll2014l l try to make an ad¬ 
vance in the treatment of the radiative feedback, mainly 
by rela xing the assumptio n of a uniform rei onisat ion field 
llMaccid et al.l 120101 : ILL De Lucia fc Helmil l2010h . These 
studies generally combine the merger-tree histories of ad¬ 
vanced N-body simulations fe.g. lSpringel et, a l. 2008) with a 
semi- analytic radiative transfer (see for instance Zahn et al.1 
dm!) describing an in-homogeneous reionisation process. 
Despite the different approaches, the population of faint 
satellites is always highly sensitive to the adopted reioni¬ 
sation model. 


We finally point out that the obser ved properties of the 
MW satellites in the Loca l Universe llSalvadori fc Ferraral 
120091 : ISalvadori et al.l 120141 ) are the result of the global in¬ 
terplay between radiative, mechanical and eventually chem¬ 
ical feedback and a more comprehensive study is necessary 
to interpret the various observational signatures left during 
galaxy evolution: the origin and properties of the observed 
population of metal-poor stars in our galaxy and the low 
metallicity tail of the MDF. 

In this paper we introduce GAMESH, a new pipeline inte¬ 
grating the latest release of cosmologic al radiative transfer 
code CRASH liGraziani. Maselli fc CiardfeORih with the semi- 
analytic model of ga laxy formatio n GAMETE, powered by a N- 
body simulation llSalvadori et al.ll2010l . hereafter SF10). By 
following the star formation, metal enrichment and photo¬ 
ionisation in a self consistent way, GAMESH is an ideal tool to 
study the effects of photo-ionisation and heating in galaxy 
formation simulations and can address many of the ques¬ 
tions discussed above. Into the bargain, the radiative trans¬ 
fer model adopted by GAMESH relies on a Monte Carlo scheme 
and not only naturally accounts for the inhomogeneities of 
the ionisation process, but also allows a deep investigation of 
other radiative transfer effects, a self consistent calculation 
of the gas temperature, and a more accurate description of 
the stellar populations responsible for the Milky Way envi¬ 
ronment reionisation. 


In the present work the GAMESH pipeline will be ap- 
plied to the Galaxy fo rmation simulation introduced in 
IScannapieco et all (l2006l l and successfully post-processed by 
the semi-analytic code GAMETE (SF10), to derive the prop¬ 
erties of the central galaxy and to compare them with the 
available observations of the Milky Way system. Although 
the mass resolution of the adopted N-body does not pro¬ 
vide a sufficient statistics of mini-halos to study the missing 
satellite problem, the robust set of results already obtained 
with the previous semi-analytic approach offers an excellent 
validation framework for the GAMESH pipeline and allows to 
appreciate the many advantages of the accurate radiative 
transfer treatment of CRASH coupled with the semi-analytic 
approach of GAMETE. A future work, already in place, will 
focus on high resolution simulations to address all the sci¬ 
entific problems discussed above with the adequate level of 
detail. 


The paper is organised as follows. In Section 2 we de¬ 
scribe the GAMESH pipeline by briefly introducing its compo¬ 
nents and by providing details on the feedback implementa¬ 
tion. The Milky Way reionisation simulation is introduced 
in Section 3 and the results discussed in Section 4. Section 
5 finally summarises the conclusions of the paper. 
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Figure 1. GAMESH pipeline logic at fixed redshift Zi. The quanti¬ 
ties Xgai, Tgai, £ g as and Tgas refer to the ionisation fractions and 
temperatures of the gas in the grid cells containing galaxies and 
in the simulated domain, respectively. The gas number density 
projected into the grid used by CRASH is indicated as n gas and 
the global set of information provided by the N-body merger tree 
as M-Tree. The quantities used by interactor Iq are the star for¬ 
mation rates (SFR) and the metallicity of the stars (Z*), while 
the computed ionisation rates and spectral shapes per galaxy and 
population are indicated as 1V 7 and ,5V. See text for more details. 


2 THE GAMESH PIPELINE 

In this section we describe the work-flow of the GAMESH 
pipeline which integrates a N-body simulation, the semi- 
analytic code GAMETE (see Section 2.1 for more details), and 
the radiative transfer code CRASH (see Section 2.2); more 
details on each pipeline module can be found in dedicated 
subsections. 

GAMESH models the galaxy formation process concate¬ 
nating a series of snapshots provided by a N-body run at 
given redshifts Zi=o,...,N- the pipeline uses the physical quan¬ 
tities calculated at Zi as initial conditions for the successive 
calculation at Zi+ 1 . Along the redshift evolution, the feed¬ 
back between star formation and RT is handled by two soft¬ 
ware modules called interactors Jo (detailed in Section 2.3) 
and 1 1 (see Section 2.4). Io transforms the galaxy star for¬ 
mation rates (SFR) predicted by GAMETE into a list of ion¬ 
ising sources for CRASH, while Ji uses the gas ionisation and 
temperature determined by the RT to establish a star forma¬ 
tion prescription for the semi-analytic model implemented 
in GAMETE. Hereafter we focus on the pipeline logic at fixed Zi 
(see Figure 1); for an easier reading we also generically refer 
to the ionisation of the gas as x SS lS , while a more specific no¬ 
tation will be adopted in Section 4 to discuss the ionisation 
fractions of hydrogen and helium. 

The initial conditions of the pipeline (ICs) are provided 
by the N-body simulation which assigns the simulation red¬ 


shift Zi to all the components, sets up the N-body merger- 
tree into 1 1 and the gas number density (n gas ) in the grid 
used by CRASH to map the physical domain. 

Once the ICs are set up, I\ starts the simulation by cre¬ 
ating a list of galaxies found in the merger-tree: each galaxy 
is identified by a unique ID and it associates the values of 
Xgai and Tg a i found in the cell of the grid containing the 
galaxy center of mass (see section 2.4 for more details). This 
list is then processed by GAMETE to establish which galaxy 
can form stars, self-consistently with the metallicity, tem¬ 
perature and ionisation of the accreting gas. 

As output of GAMETE we obtain a sub-sample of star 
forming galaxies, their SFR,, stellar metallicity and popula¬ 
tion type. Io converts this sample into a list of CRASH sources 
by evaluating the galaxy positions on the grid, their spec¬ 
trum integrated ionisation rate 1V 7 and the spectral shape 
S v . A stellar synthesis database has been implemented in 
Io (see Section 2.3) to derive A 7 and S v from the stellar 
metallicity and the SFR. 

Once the properties of the radiating galaxies are estab¬ 
lished, the radiative transfer simulation starts propagating 
photons for a simulation duration corresponding to the Hub¬ 
ble time separating two snapshots, and it obtains the gas 
ionisation a; gas and temperature T gas at redshift Zi. These 
quantities are finally used for the subsequent redshift Zi+ 1 , 
by repeating the same algorithm. 

Before moving to a more detailed description of the 
GAMESH components, we want to emphasise the advantages of 
our approach. We first point out that the use of CRASH allows 
us to follow with accuracy the process of reionisation of the 
MW progenitors along the redshift evolution, by taking into 
account the intrinsic inhomogeneities due to the gas clumps 
and by calculating the temperature history self-consistently. 
Moreover, radiative transfer effects, which generally lead to 
spectral hardening that may preferentially heat gas in over- 
dense regions, will be accounted for by the RT simulation 
without pre-assuming any propagation model. 

The adoption of GAMETE allows us to fully explore the 
interplay between reionisation, radiative feedback and chem¬ 
ical evolution which can have very interesting and testable 
consequences (Schneider et al. 2008). Using this approach, 
we can break the degeneracy between having gradual sup¬ 
pression and constant star formation efficiency or sharp sup¬ 
pression and mass-dependent star formation efficiency. 

In the following sections more details on the pipeline 
components are provided. 


2.1 GAMETE 

GAMETE (jSalvadori. Schneider fe Ferraral 12007 ) is a data- 
constrained semi-analytic model for t he formatio n of the 
Milky Way and its dwarf satellites (ISalvadori fe Ferraral 
1200911 . The algorithm has been designed to study the prop¬ 
erties of the first stars and the early chemical enrichment 
of the Galaxy. The evolution of gas and stars inside each 
galactic halo of the merger-tree hierarchy (group of Dark 
Matter (DM) particles) is traced by assuming the following 
hypotheses: 

• the gas is assumed of primordial composition at the 
highest redshift of the merger tree; 
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• in each galaxy the SFR is proportional to the cold gas 
mass; 

• the contribution of radiative feedback is accounted for 
by adopting different prescriptions, depending on the prob¬ 
lem at hand. When the code runs in stand-alone mode (i.e. 
not in pipeline) GAMETE assumes instant reionisation (IReion, 
SF10), i.e. stars form only in galaxies of mass Mh > M±(z) = 
3 x lO 8 M 0 (l + 2 )- 3 / 2 (Mh > M 3 o(z) = 2.89 x Mi(z)) prior 
to (after) reionisation (assumed complete at z = 6). When 
GAMETE runs coupled with CRASH, star formation is not regu¬ 
lated by the halo mass but by exact radiative feedback effects 
as predicted by the RT simulation and detailed in Section 
2.4. In this case the star formation efficiency in mini-halos 
(i.e. halos with T v j r < 2 x 10 4 K) is decrease d oc T)” 3 as al¬ 
ready detailed in lSalvadori fc Ferrarai (2012), to mimic the 
effects of a LW background. 

Although recent versions of GAMETE (Salvador i et al.l 
l20l4 fde Bennassuti et al1l2014} ) allow to follow the Pop III 
/ Pop II transition by accounting for detailed chemical feed¬ 
back (including the dust evolution) and exploring different 
Pop III IMFs and stellar lifetimes, as well as the effects of 
the in-homogeneous distribution of metals, the first release 
of GAMESH implements the simplified version of the chemical 
network as detailed below, for a better comparison with the 
results in SF10. 


• As assumed in the so cal led ’critical metallicity scenario’ 
llSchneider et alJl2002j . l200fj ) low-mass stars form following 
a standard, Salpeter-like Initial Mass Function (IMF), when 
the metallicity of the gas is Z > Z cr = 10” 4 Zq; when Z < 
Z cr we assume that Pop III stars form with a mass m Pop m = 

200Mq . 

• The enrichment of gas within the galaxies and in the 
diffused MW environment (hereafter Milky Way Environ¬ 
ment MWenv) is calculated by includin g a si mple descriptio n 
of sup ernova (SN) feedback llSalvadori. Ferrara fe Schneider! 
l2008h . By adopt ing the Instantaneous Recycling Approxi¬ 
mation ljTinslevl 119801 ) we assume that the gas is instanta¬ 
neously and homogeneously m ixed with the atomic metals 


Salvadori. Schneider fe Ferr ara ( 2003) ; ISalvadori et al.1 
20icl ) for a discussion on the implications of this choice). 


The interplay between N-body scheme and semi- 
analytic calculation, along the redshift evolution of the sim¬ 
ulation, is detailed below. At each time-step we distribute 
the mass of gas, metals and stellar component of each halo 
among all the DM particles; this is used as ICs of the suc¬ 
cessive integration step. Metals ejected into the MWenv are 
treated in the same way so that the newly virialising halos 
have a chemical composition which depends on the level of 
enrichment the environment in which they are formed. To 
recover the spatial distribution of the long-living metal-poor 
stars at z = 0, we also store their properties per DM particle. 


2.2 CRASH 


CRASH (|Ciardi et a 

. 200ll; Maselli, Ferrara & Ciardi 

2003; 

Maselli & Ferrara 

2005!: Maselli. Ciardi & Kanekar 

200£ 

Graziani, Maselli & Ciardi 20131) is a Monte Carlo 

sased 


scheme implementing 3D ray tracing of ionising radiation 
through a gas composed by H and He and atomic metals 
(e.g. C, O, Si). The code propagates ionising packets with 


N~f photons per frequency along rays crossing an arbitrary 
gas distribution mapped onto a Cartesian grid of TV 3 cells. 
At each cell crossing, CRASH evaluates the optical depth r 
along the casted path and the amount of absorbed photons 
N a hs = IV 7 ( 1 — e T ). Nabs is used to calculate the ionisation 
fractions of the species *g as = x± G (ithii, £Heii, *Hein) and 
the gas temperature Tgas of the crossed cell. 

The initial conditions of a CRASH simulation are assigned 
on a three-dimensional Cartesian grid and a cosmological 
box of linear size Lb, as detailed in the following list: 

• the number density of H (hh) and He (riHe), the tem¬ 
perature of the gas T gas and the gas ionisation fractions x± 
at the initial simulation time tcQj 

• the number of point sources emitting ionising radiation 
(N a ), their position in the Cartesian grid, their ionisation 
rate (N y in photons s” 1 ) and finally their spectral energy 
distribution (SED, £„ in erg s” 1 Hz” 1 ). The SED is assigned 
as an array of frequency bins providing the intensity of the 
radiation; 

• the duration of the simulation ( tf ) and a pre-assigned 
set of simulation times tj G {to,...,i/} used to store the 
relevant physical quantities; 

• the intensity and SED of a UVB, if present. 

The interested reader can find more technical informa¬ 
tion about the lat est i mpl ementatio n of CRASH and references 
on its variants in lGraziani, Maselli fe Ciardil ll2013l) . 


2.3 Radiation sources 

As explained above, CRASH requires the list of all the ion¬ 
ising sources present in the box grid. This is provided in 
the GAMESH framework by the GAMETE-to-CRASH module Iq, 
responsible of converting the properties of the star forming 
galaxies (Pop II/Pop III SFR, Z») into CRASH sources with 
specific spectral properties. 

Jo first maps the center of mass of each galaxy onto grid 
coordinates and then converts the star formation rates into 
spectrum-integrated quantities iV 7 depending on the stellar 
population. 

For Pop II stars we c alcu l ate ionisatio n rate s and spec- 
tral shapes accordingly to lBruzual A. fc Chariot! j| 1993 ! and 
we assume a IMF in the mass range [0.1 — 100 ]Mq. A 
different spectral shape and ionisation rate is then associ¬ 
ated with each of the Pop II star forming galaxy depend¬ 
ing on its population lifetimes t* and stellar metallicity 
Z *. The spectrum and ionisation rate is derived from a 
grid of pre-computed spectra integrated in specific lifetime 
bins t» G {0.001, 0.01, 0.1, 0.4,1.0, 4.0,13.0} Gyr and stellar 
metallicity Z* G {0.005, 0.2, 0.4,1.0, 2.5}Z 0 . 

For Pop III stars we assume an ionisation rate per sol ar 
mass, N y = 1.312 ■ 10 48 [photons s _1 /Mg] llSchaereij 2002 :• 
corresponding to a stellar mass M* ~ 2OOM 0 , averaged on a 
lifetime of about 2.2 Myr. The SED of Pop III stars is simply 
assumed as black body spectrum at temperature Tbb = 
10 5 K. 


1 fi'gas is calculated by projecting the particle distribution of the 
simulation on a Cartesian grid of N% cells and by deriving the 
gas component via the universal baryon fraction. 
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2.4 Radiative feedback 

The radiative feedback on the star formation is implemented 
by the pipeline module I\. This module is responsible for set¬ 
ting up the ionisation and temperature (a; ga i, ^gal) of the gas 
in the galaxies, successively processed by GAMETE to establish 
their star formation. 

At the first redshift zo it is simply assumed that the 
medium is fully neutral (i.e. :r ga i = 0) and T ga i = To(l + 
zo), where To is the value of the CMB temperature at z = 
0. During the redshift iterations Ji computes a; ga i, ^gai by 
first finding the cell containing the galaxy center of mass 
in the CRASH grid. As second step, it evaluates which of the 
surrounding cells best describe the environment from which 
the galaxy can get the gas to fuel star formation. 

The effective environment scale from which cold gas is 
fuelled into star forming galaxies is crucial to calculate the 
values of * ga i, T ga i and to apply the radiative feedback. As 
in our pipeline the gas distribution surrounding halos re¬ 
lies on the spatial resolution of the RT (i.e. the cell size in 
CRASH: AT = L b /N c ), we compare the virial radius R v i r 
of the galactic halo with AL/2. We simply assume that if 
2R v i r /AL ^ 0.1, the galactic environment is mainly set 
up in the cell containing its center of mass and we as¬ 
sign a-gai = a; C eii, T ga i = T C eii. When, on the other hand, 
2R v i r /AL >0.1 the gas reservoir of the galaxy could ex¬ 
tend to the surrounding cells and then T gai is assigned to 
their volume averaged value; the value of x ga i remains in¬ 
stead the one taken from the central cell. It should be noted 
though, that the threshold value 0.1 depends on the reso¬ 
lution of the CRASH grid and must be tuned for increased 
resolutions grids, when necessary 0. 

Once the temperature and ionisation fractions affecting 
the galactic star formation have been assigned, we first use 
* ga i to evaluate the mean molecular weight fi of the gas and 
then to calculate the viri al temperature of the g alactic halo 
T v i r (see the formulas in lBarkana fc Loebll200ll l. 

Star formation in the galaxy is finally admitted in 
GAMETE if the condition T ga i < T v i r applies, allowing accre¬ 
tion of cold gas in the galaxy to form stars. In the next ver¬ 
sion of GAMESH we will refine the star formation prescription 
also accounting for gas cooling, metallicity dependence and 
dynamical time scales following the approach described in 
iNoh fe McQuinnl (j2014h : also see the Introduction for more 
references and details. 


3 MILKY WAY REIONISATION SIMULATION 

Here we describe the set-up of our Milky Way reionisation 
simulation. After a brief description of the N-body simu¬ 
lation and its parameters we detail the radiative transfer 
assumptions a nd the initial condi t ions. 

Following IScannapieco et ahl (j2006l ). we adopt a cold 
dark matter (ACDM) cosmological model with h = 0.71, 


2 Improving the accuracy of the galactic environment description 
would ideally require to selectively resolve the halo structures (i.e. 
a cell size AT << R v i r ), as obtained by using cell refinements 
around the galaxies. The spatial resolution adopted in the paper 
simulation is documented in Section 3.2 
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Qoh 2 = 0.135, Da = 1 - D 0 , tl b h 2 = 0.0224, n = 1 and 
as = 0.9. 


3.1 The N-body simulation 


To study the Milky Way reionisation we adopt the 
snapshots and merger tree from a cosmological N-body 
simulation of the Milky Way-sized galaxy halo carried 
with GCD+ dKaw ata fe Gibsonl 200.3al ) . which were used in 
IScannapieco et a l.l (20o dl. The detai l of the simulation is 
described in Scannapieco et ah] (l2006l l: here we briefly sum¬ 
marise the simulation parameters necessary to understand 
the global GAMESH run. 

The mass adopted for the N-body particles is M p = 
7.8 • 10 s Mg and the softening length used is 540 pc. The 
simulated system consists of about 10 6 particles within the 
virial radius R v i r = 239 kpc for a total virial mass M v i r = 
7.7 ■ 10 11 Mq. Note that the virial mass and virial radius 
estimated by observations for the Milky W ay are M v j r = 
10 12 Mq and R v i r = 258 kpc respectively llBattaglia et al.l 
l2005lf . 

_ Using a multi-resolution technique llKawata fe Gibsonl 

2003b), the initial conditions at z = 56 are set up to resolve 
the region in a sphere within the radius of R = 4R v i r with 
the high-resolution N-body particles of mass M p = 7.8 • 
1O 5 M0 and with the softening length of 540 pc; the other 
region is resolved instead with lower resolution particles. 

The simulation snapshots are provided with regular in¬ 
tervals every 22 Myr between z = 8 and 17 and every 110 
Myr for z < 8. The virialised DM halos are then identified 
by friend-of-friend group finder by assuming a linking pa¬ 
rameter b = 0.15 and a threshold number of particles of 50; 
the resulting minimum mass halo is Mh = 3.75 x 10 7 Mq. 

A validating N-body simulation at lower resolution, also 
including a star formation recipe, has also been performed 
to confirm that t he assu me d in itial condition lead to a disc 
formation fsee lBrook et al.ll2007l l. 


3.2 The RT set-up 

The radiative transfer simulation is performed at each red¬ 
shift Zj, on a box size of 2/i _1 Mpc comoving, by emitting 
N p = 10 6 photon packets from each star-forming galaxjtj. 
The box is mapped on a Cartesian grid of N c = 128 cells per 
cube side, providing a cell resolution of AT ~ 15.6 h~ x kpc. 
The gas in the box is assumed of cosmological composition 
(92% H and 8% He) and we neglect the metal ions and the 
effects of the metal cooling around the star-forming halos. 
We defer to a future application the extension of GAMESH 
to the metal ions. For a gas of cosmological composition 
(n H i ~ 92% n gas ) the ionising spectrum is usually assigned 
in the energy range 13.6 eVsJ _E 7 ^ 140 eV. Note that the 
spectral database could be extended at any time to include 
harder spectra accounting for sources emitting higher en¬ 
ergy photons. Note also that the CRASH calculation adopts 
the same cosmological parameters as in the N-body simula¬ 
tion to maintain consistency across the pipeline life-cycle. 


3 Note that according to the pipeline algorithm each source could 
have a ’a priori’ different spectrum assigned by 1$ as function of 
its lifetime and stellar metallicity. 
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z 

Figure 2. Redshift evolution of the volume averaged ionisation 
fractions and temperature of the gas, down to z ~ 5. From top 
to bottom the values of mii, *Heii, aiHeiii, T[10 4 ] K are shown in 
solid green, blue, black and red lines respectively. 


The convergence of the Monte Carlo sampling is guar¬ 
anteed by an identical run with a factor of ten larger N p 
and providing the same numbers within the fourth decimal 
of the volume average ionisation and temperature created in 
the box. 

The simulation starts at zo ~ 1(0 with uniform tem¬ 
perature To ~ 46 K and we assume that the reionisation is 
completed when the volume-averaged ionisation fraction of 
the hydrogen reaches the value xmi ^ 0.995 against the gas 
recombination. Due to the large excursion in Hubble time, 
a gas recombination CASE-B is adopted in the simulation 
and the diffuse re-emission is neglected. 

In the absence of a model for the UV background on 
the scale of the MW formation, we avoid using an exter¬ 
nal cosmological UVB as done in other works present in the 
literature, and prefer to mimic a flux entering our box by 
applying periodic boundary conditions to the escaping ra¬ 
diation. While this choice guarantees that the photon mean 
free path is conserved when the box ionisation is advanced 
and the gas becomes transparent to the hydrogen-ionising 
radiation, it should be noted that we do not account for an 
external ionising flux likely emitted by the background of 
QSO sources established below z ~ 4. 


4 RESULTS 

In this section we discuss the results of our simulation, also 
comparing them with trends obtained adopting the assump¬ 
tion of instant reionisation at 2 ~ 6 as discussed in SF10. 
Hereafter, for an easier reading, we omit the label "gas" in 
the variable names referring to ionisation fractions and tem¬ 
perature of the gas and we comment on the single ionised 
species x± G (ihii, ZHeii, ZHein). 

4.1 Reionisation and temperature histories 

Here we investigate the evolution in redshift of the volume 
averaged ionisation fractions and temperature (in this sec¬ 
tion Xi(z),T(z)) resulting from the radiative transfer simu¬ 
lation. 

In Figure 2 we show the redshift evolution of iehii (top 
panel) and * H eii, *Hein (second and third panels from top) 
down to 2 = .0. Xhii remains below Xyhi ~ 0.1 when 2 > 12, 
afterwards it rapidly rises up to *hii ~ 0.5 at 2 ~ 11 
and reaches x H n ~ 0.9 when z ~ 10. The sudden increase 
in the hydrogen ionisation fraction in the redshift range 
10 < 2 ^ 12 can be ascribed both to the rising in the co¬ 
moving SFR (Figure 4) and in the total number of emitting 
galaxies (Figure 5) but more likely to an advanced stage of 
overlapping of ionised regions, which makes the reionisation 
process intrinsically non-linear. 

Below 2 ~ 10 the hydrogen reionisation wades through 
and *hii increases up to 0.99 by 2 ~ 9 and up to 0.999 by 
2^5. Even if not visible in the Figure, we point out that 
below 2 ~ 8.5 the value of the third decimal of i H ii continues 
to oscillate in 0.997 < Xnn ^ 0.999 due to the continuous 
unbalance in the gas ionisation, created by the ongoing gas 
collapse that enhances the recombination rate at the cen¬ 
ter, and by the in-homogeneous reionisation/recombination 
occurring in the volume. We consider the hydrogen reioni¬ 
sation completed at z ~ 6.4 when Xau > 0.995, neglecting 
these fluctuations. 

The second and third panels from top, show the ion¬ 
isation histories of helium. It is immediately evident that 
the He i reionisation can be sustained only by the Pop II 
stars dominating below 2 ~ 12 while the reionisation of 
He ii never occurs at the considered Milky Way scale (here 
2/i -1 Mpc comoving) without the contribution of an exter¬ 
nal background of harder spectra (e.g. QSO). In fact, while 
the redshift trend of a^eii follows the one of the hydrogen, 
a?Hem never reaches values higher than ^Hein = 0.002 along 
the entire redshift range. As also pointed out for Zhii, the 
volume averaged ionisation fraction of He ii continues to os¬ 
cillate in 0.97 < ZHeii < 0.99 below 2 ~ 8.5, confirming 
that a harder spectral component is required to sustain full 
helium ionisation against gas recombination. Note that a 
bump in the value of aiHein (third panel from top) occurs 
in the redshift range 12 < 2 < 15, tracing the presence of 
Pop III stars (see also Figure 4 ) with their harder spectral 
shapes. Note also that the volume averaged value of a^em 
remains always negligible over the entire simulation: after 


5 Below this redshift a uniform ionised pattern of hydrogen is 
4 This is the first redshift in which the applied FoF is able to created and the values are less indicative. Note that we cannot 

identify the first halos. This redshift value depends on mass res- discuss the helium reionisation due to the absence of a harder 

olution of the simulation. radiation background. 
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Figure 3. Slice cuts of n ga s (first column), ojhii (second column), T (third column) at three different redshifts: z ~ 12 (top panels), 
z ^ 11 (middle panels), z ~ 6 (bottom panels). The value of the fields in the plane are represented as colour palettes and the distance 
units are expressed in number of grid cells along the x and y axis (~ 15.6/i —1 kpc per cell). In the first-column panels star forming 
galaxies are symbolised as blue crosses when forming Pop III stars, while as gold-yellow asterisks for Pop II stars; the sizes of the symbols 
are scaled with the galaxy ionisation rate TV 7 
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the transition from Pop III to Pop II stars the presence of 
full ionised helium is so closely confined to the sources that 
is not visible in the panel, except at z < 7, when hydrogen 
ionisation is complete and harder photons (f? 7 > 54.4 eV) 
are free to cross the box many times and to create a very 
small ionisation fraction of fully ionised helium. 

The redshift evolution of the volume-average tempera¬ 
ture is shown in the bottom panel of Figure 2 in [10 4 K] scale. 
Before 2 ~ 10, the temperature rises from the initial value 
T ~ 50 K up to T ~ 10 4 K; note that the rapid evolution of 
Pop III stars in this redshift range is so small (not visible in 
the plot) that it cannot change the global trend of the aver¬ 
age temperature. In the redshift interval 9 < 2 < 10 the gas 
stabilises at a average photo-ionisation equilibrium tempera¬ 
ture ofT~ 10 4 K, before its average value starts decreasing 
by some 10 3 K during the successive evolution. Oscillations 
of the order of 10 3 K are present below z ~ 5 (not shown 
in the plot), due to the competitive effects of gas clumping 
toward the center, adiabatic cooling due to the cosmologi¬ 
cal expansion, and finally in-homogeneous distribution and 
properties of the ionising galaxies. Even if not visible in this 
plot we report a final average temperature of T = 7.9 x 10 3 
K at z — 0 but we remind that this value is certainly un¬ 
derestimated in our simulation due to the absence of high 
energy photons from the QSO background established on 
the large scale below a ~ 4. 

To visually illustrate how the reionisation process 
evolves in space and redshift, we show in Figure 3 three slice 
cuts of our volume (of 2 hr 1 comoving Mpc side length) taken 
at redshifts z ~ 12,11,6 (see panels from top to bottom). 
Panels in the first column show log(n gas (z)) as gray gradi¬ 
ent from black ( n ga3 ~ 2 x 10 -5 cm -3 ) to white ( n gas ~ 0.4 
cm' 3 ); in the second column Xmi(z) is shown by a colour 
palette from violet (ihii = 10 -5 ) to orange (xmi = 1.0), 
and finally the pattern of T(z) is drawn in the third column 
with a different colour palette from white (T ~ 50 K) to 
dark red (T ~ 1.9 • 10 4 K). In all the panels, the distances 
are shown in cell units (~ 15.6/i _1 kpc per cell) along the x 
and y directions. 

Superposed to the number density maps we also show 
the star forming galaxies (SFG) present in the surface of 
choice: galaxies forming Pop III stars are symbolised as blue 
crosses (top panel), while gold-yellow asterisks are associ¬ 
ated with the presence of Pop II stars. Also note that the 
size of each symbol is scaled with its IV 7 and the position of 
star forming galaxies perfectly correlates with the gas fila¬ 
ments, falling toward the center where the central galaxy is 
progressively forming along the redshift evolution. 

As consequence of the central position of the ionising 
sources, the reionisation proceeds inside-out. At high red¬ 
shift (top panels), the low density regions surrounding the 
center of the box (black regions, first column) are easily 
ionised by the sources (orange pattern, second column). At 
the center of the image note that high-density clumps (white 
regions) are able to suppress photo-ionisation with a high 
gas recombination rate (magenta patterns). The average 
gas temperature (third column) tends to stabilise around 
T ~ 1.4 x 10 4 K in the regions involved by the overlap 
of ionised bubbles (yellow and red), mainly driven by the 
contribution of Pop III stars. T rapidly decreases down to 
T ~ 10 2 K at the border of the ionised region (blue pat¬ 
terns) not reached by stellar radiation. Consistently with 


the ionised pattern, few high-density structures remain self- 
shielded and settle at very low temperature at the center 
image. Note that the presence of helium in our simulation 
significantly changes the gas photo-heating enhancing the 
contrasts between over-dense regions and voids. At z ~ 11 
(middle panels) the reionisation process is already quite ad¬ 
vanced and the ionising radiation involves a large fraction of 
the box, including the under-dense regions far away from the 
center (green and yellow patters). Voids are rapidly ionised 
and their temperature stabilises around T ~ 6 • 10 3 K, well 
below the central value T ~ 10 4 K. Finally, at 2 ~ 6 the en¬ 
tire volume reaches photo-ionisation equilibrium and its av¬ 
erage temperature stabilises around T ~ 8.0 • 10 3 K but the 
spatial pattern continues to show large variations from the 
gas filaments (branching off the central MW-type galaxy) 
and from the under-dense regions. 

Note that at the final redshift (z ~ 6.1) the central 
galaxy is seen as single source with highest ionisation rate 
of about iV 7 = 2 • 10 54 photons/s and SFR ~ 13Mg/yr 
(the total SFR of the emitting objects is ~ 20Mq/yt ). In 
the other galaxies laying on the plane, the star formation is 
suppressed by radiative feedback (see Section 4.2.1 for more 
details) 

As final comment, we point out that our simulation 
lacks the ionising contribution of an external UV back¬ 
ground, especially at the helium ionising frequencies. While 
periodic boundary conditions are adopted to preserve the 
photon mean free path at all frequencies, the external UVB 
may play a relevant role because of the small size of the 
simulation box and its coarse mass resolution. On the other 
hand, the inclusion of a UVB is neither trivial nor straight- 
forward. First, the ava ilable UVB models (see for example 
lHaardt fe Madaull2012l and references therein) are computed 
and calibrated on the large scale structure. Therefore, on the 
scales of MW formation, both the intensity and the spectral 
shape of the UVB may be modified by RT effects. By means 
of the radiation tracking features of CRASH3, we will investi¬ 
gate these effects in a future study. Here we only point out 
that an external UVB could only impact hydrogen ionisation 
above z ~ 11. In fact, Fig. 2 shows that below this redshift 
the internal flux is sufficient to sustain hydrogen ionisation. 

4.2 Feedback on star formation 

The redshift evolution of the comoving SFR is shown in Fig¬ 
ure 4, where the red lines refer to the contribution from Pop 
II stars (highest curve) and the blue lines to Pop III. Dotted 
lines, with same colours, show the values obtained by adopt¬ 
ing the instantaneous reionisation prescription (See Section 
2.1). In the redshift interval 12.4 < z < 16 the total SFR of 
our simulation (solid lines) increases up to 1.7Mq /yr thanks 
to both Pop III and Pop II stars, but note that at the highest 
redshift z ~ 16.5 (when the gas is mostly primordial) only 
population III objects contribute to the total SFR. Feedback 
from Pop III forming galaxies, both radiative and chemical, 
acts efficiently on the surrounding medium: their hard spec¬ 
tra locally increase the MWenv temperature inhibiting star 
formation in unpolluted clouds, while mechanical feedback 
rapidly pollutes the MWenv up to the critical metallicity 
Ecr — 1O” 4 Z 0 . Both effects contribute to their sudden dis¬ 
appearance, triggering the formation of Pop II stars below 
z ~ 12.4. At lower redshift Pop II stars cause an increase of 
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Figure 4. Star formation rate for Pop II (solid red lines) and Pop 
III (solid blue) as function of z. As reference the values obtained 
with the instant reionisation approximation (IReion) are shown 
dotted lines and same colours for the two populations. The value 
of the SFR estimated by observation at z = 0 is reported as error 
bar. See text for more details. 

the total SFR up to a value of SFR ~ 25Mq/yi at the peak 
redshift z ~ 4. The successive evolution proceeds with an 
irregular but progressive decrease down to a value of SFR 
~ 1.8 Mp)/yr at z = 0, in good agreement with the observed 
data Ichomiuk fe Povichll20li] ). 

It is interesting to compare these results with the ones 
obtained with the instant reionisation (IReion) approxima¬ 
tion for both populations. We remind that in IReion approx¬ 
imation GAMETE accounts for feedback by assuming that only 
Lyman-a halos having Mh > M± can form stars at z > 6, 
while at z ~ 6 the volume is suddenly ionised and a different 
threshold applies. For z < 6 only halos with circular veloc¬ 
ity greater than 30 kms -1 (i.e. Mh > M 30 ) can trigger the 
star-formation (see SF10). 

By comparing the solid and dotted lines in Figure 4 we 
immediately conclude that the suppression of star formation 
operated by the IReion is quite appropriated to reproduce 
the SFR of Pop II stars in a wide redshift range but it over¬ 
estimates the SFR at z = 0 by 15%. Note also that the 
stronger radiative feedback effects at high redshift in Ireion 
lower the SFR of Pop III stars and delay the transition be¬ 
tween the two stellar populations. 

4-2.1 Feedback statistics 

The effects of the radiative feedback on the statistics of star 
forming galaxies are shown in Figure 5 as function of redshift 
z. The solid blue line indicates the number of SFG computed 
in our simulation, while the dotted red line refers to the 
IReion case. The total number of DM halos found in the 
N-body simulation in use is shown as dashed black line. 
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Figure 5. Number of star forming galaxies as function of z in 
the simulation (solid blue line). The dotted red line refers to the 
IReion case, while the dashed black line shows, as reference, the 
total number of DM halos found in the N-body simulation. 

By comparing solid and dotted lines, we first note that 
at the current resolution of the simulation, which does not 
resolve the full range of masses that sample the mini-halos 
regime, the condition Mh > M 4 could appear quite appro¬ 
priate to reproduce the statistics of the galaxies with sup¬ 
pressed star formation at very high redshifts (z > 10). On 
the other hand, a careful comparison with Figure 4 immedi¬ 
ately shows a different star formation rates of Pop III galax¬ 
ies (compare blue lines in Figure 4). 

This discrepancy can be ascribed to a combination of 
reasons which involve radiative, mechanical, and chemical 
feedback. Also note that the IReion model cannot account 
for in-homogeneous feedback on high-redshift low mass ob¬ 
jects (which are likely to host Pop III stars) but indiscrimi¬ 
nately applies to all the galaxies with Mh > M 4 present in 
the volume. At high redshift this assumption is particularly 
incorrect because the SFG are still quite separated in space 
and then their radiative feedback acts very selectively, only 
affecting their surrounding neighbourhoods. See for exam¬ 
ple how dis-homogeneous is the bubble overlap in the third 
panels of Figure 3 and how many regions of the simulated 
volume remain unaffected by photo-heating. 

Note that the current mass resolution of the N-body 
simulation does not allow to draw robust conclusions on the 
radiative feedback effects at high-redshift. This will be the 
subject of a future dedicated study (de Bennassuti et al. 
2015, in prep.). 

At lower redshifts the efficiency of the radiative feed¬ 
back increases with the progress of reionisation in both mod¬ 
els. In the redshift interval 6 < z < 9 the GAMESH run pre¬ 
dicts that around 50% of the galaxies have been affected 
by radiative feedback and this value increases up to 65% 
by z ~ 4, with a jagged trend due to both oscillations in 
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the number of potential galaxies (dashed black line) and the 
in-homogeneous heating in space induced by the RT0. Be¬ 
low 2 = 4 the suppression of star formation in small galaxies 
progresses rapidly allowing few remaining objects in the box 
to form stars, as observed in the Milky Way environment. 

Note the difference in the statistics of SFG below 2 ~ 6 
obtained in the IReion approximation. In fact a realistic ra¬ 
diative feedback operates with a more gradual suppression 
in time, due to the in-homogeneous nature of the bubble 
overlap in space. In 3 < 2 < 6 this difference induces an 
error in the number of SFG as high as 40% between the 
two cases (compare solid blue line with dotted red line). On 
the other hand this statistical discrepancy between models 
is not reflected in both the total SFR (Figure 4) and the 
MWenv enrichment (see Figure 8) which are very similar 
in this redshift range, also indicating that the high-mass 
galaxies surviving the IReion prescription provide a domi¬ 
nant contribution to these quantities. 

Figure 6 visually shows how the in-homogeneous radia¬ 
tive feedback operates in space, by comparing the same slice 
cuts and redshift sequence commented in Figure 3. Here, on 
top of the gas number density field we have traced few solid 
line iso-contours of the gas temperature: T ~ 100 K in cyan, 
T ~ 4x 10 3 K in blue, T ~ 10 4 K in green, T ~ 1.3xl0 4 K as 
orange and finally T ~ 1.5 x 10 4 K in dark-red. All the star 
forming galaxies are symbolised by yellow asterisks, while 
the galaxies in which star formation has been suppressed by 
radiative feedback are shown as black filled circles. 

By comparing the three panels it is immediately evident 
how the inhomogeneities in the reionisation process induce 
large differences in randomly suppressing star formation (left 
and middle panel) before a UVB is fully established (right 
panel). At high redshifts few suppressed SFGs (black cir¬ 
cles) are trapped in regions with T ~ 10 4 K, not necessar¬ 
ily connected with the emitting galaxies. Note for example 
the galaxy on the right side of the middle panel: even far 
from the center it is surrounded by gas with high temper¬ 
ature, that likely originates from more complicated three- 
dimensional effects created by the RT. On the other hand 
at 2 « 6 (right panel) the medium is already pervaded by a 
quasi-homogeneous UVB and then all the galaxies below a 
critical mass are easily suppressed in the entire domain. To 
study if a critical mass for the suppression of SFGs is statis¬ 
tically restored against the non-linearity induced by the RT, 
we should rely on a higher resolution simulation providing a 
better statistics for both Lyman-a and H 2 -cooling halos, as 
well as better spatial resolution. With this study we will able 
to provide a numerically motivated recipe for semi-analytic 
models. We then defer this point to the next work (de Ben- 
nassuti et al. 2015, in prep.) and limit the discussion to this 
figure, as an illustrative example of how radiative feedback 
implemented in GAMESH works. 


6 As commented in the previous paragraph, only with a higher 
mass and spatial resolution simulation we can fully understand 
the impact of an in-homogeneous reionisation process on the 
statistics of SF galaxies 
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Figure 7. Percentage of star forming galaxies as function of z 
(solid black line) and percentage of galaxies affected by radiative 
feedback in their own cell (dashed red line). The blue dotted line 
shows the percentage of galaxies forming mainly Pop II stars. 

4-2.2 Local versus environmental feedback 

As introduced in Section 2.4, the radiative feedback defined 
in GAMESH depends on the spatial extension of the galactic 
environment feeding cold gas for star formation. Depending 
on the virial radius of the dark matter halo hosting each 
galaxy and the spatial resolution of the RT grid, the feedback 
could act locally to the sources (i.e. accounting just for the 
temperature in the galaxy cell (2R v i r ^ 0.1AL)) or it could 
act globally involving larger scales (i.e. cells surrounding the 
one containing the galaxy) when 2 R v i r > 0.1 AL. 

The statistics of these two mechanisms can help to un¬ 
derstand the importance of ’local’ versus ’global’ feedback 
along the progress of reionisation, establishing a H-ionising 
uniform UVB. In Figure 7 we show the percentage of SFG 
(solid black line) together with the percentage of galaxies 
that are influenced by ’local’ feedback (dashed red line) Q. 
To understand which stellar population is affected by local 
or environmental effects we also show in blue dotted line the 
percentage of galaxies forming Pop II stars. 

As already commented in the previous figure, the per¬ 
centage of SFG (solid black line) rapidly decreases from high 
to low redshifts following the progress of reionisation with 
a spiky and irregular trend down to 2 ~ 4. The presence of 
a "plateau" in the redshift interval 6 < 2 < 9 could mark 
the existence of an extended reionisation epoch in which the 
gas is kept in photo-ionisation equilibrium in the entire box 
at T ~ 10 4 K and the hydrogen ionisation fraction increases 
from 0.99 up to 0.999. The dashed line indicates that iso¬ 
lated, high redshift galaxies (2 > 12.4) are quite insensitive 

' These number are calculated with respect to the total number 
of halos found in the simulation. 
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Figure 6. Slice cuts of the gas number density n ga s with superposed temperature contour-plots: T ~ 100 K (cyan), T ~ 4 x 10 3 K 
(blue), T rsj 10 4 K (green), T ~ 1.3 X 10 4 K (orange) and T ~ 1.5 x 10 4 K in dark-red. All the star forming halos in the plane are 
represented by yellow asterisks, while black dots indicate halos in the plane in which star formation is suppressed by radiative feedback. 


to their environment and the radiative feedback plays only 
a local role. A quick comparison with the blue dotted line 
clearly shows that these galaxies mainly form Pop III stars. 
After the population transition, Pop II forming galaxies are 
instead progressively affected by their environment: below 
2 ^ 8 , both the increase of the virial radii of the sources and 
a more uniform ionising background make the environment 
progressively more important and the feedback acts glob¬ 
ally. At these redshifts the gas temperature is sustained by: 
i) the few satellite galaxies surviving radiative feedback, ii) 
the global UVB entering the box, and iii) the predominant 
emission from the central MW-type Galaxy. It is in fact clear 
from Figure 7 that less than 40% of galaxies is capable to 
accrete gas from the surrounding reservoir to fuel the star 
formation process below 2 ~ 4. 

4.3 Interplay with chemical feedback 

In this section we briefly investigate how the reionisation 
process can leave specific signatures in the chemical evolu¬ 
tion of the Milky Way environment and in the final prop¬ 
erties of the MW at 2 = 0 and, in particular, of the most 
metal-poor stars observed in the Galactic halo. 

In Figure 8 we show the redshift evolution of the MWenv 
metallicity ZMWemlE predicted by GAMESH (solid black line) 
and the corresponding values obtained by assuming instant 
reionisation (dotted red line). As already noted comment¬ 
ing Figure 4 and 5, the star-formation history of Pop III 
stars is remarkably different when an accurate model for 
the radiative feedback is taken into account and their rapid 
burst, predicted by our simulation, is clearly imprinted in the 
metallicity of the diffuse, external medium. In fact, ZMWenv 
experiences a rapid increase at high redshift, reaching the 
critical value Z cr = 10 ~ 4 Zq already at 2 ~ 14. 

8 This value is defined as the metallicity of the medium surround¬ 
ing the collapsed halos and it is calculated as ratio of the mass of 
metals over the total mass of the gas. 



Z 


Figure 8. Evolution in redshift of the Milky Way environment 
metallicity ZMWenv (solid black line) in solar metallicity units. 
The values obtained by assuming instant reionisation are shown 
in dotted red line. The green dashed line shows, as reference, the 
value adopted for the critical metallicity to form Pop II stars. 

After a flat evolution corresponding to the Pop III to 
Pop II transition the increasing trend is restored by the in¬ 
crease in SFR of Pop II stars. As for the SFR (see Figure 4), 
ZMWenv becomes consistent with the value predicted by the 
IReion case in the redshift interval 3 < 2 ^ 8 after which 
it flattens. Note that a 15% difference in the SFR at 2 = 0 
due to an extended reionisation is reflected in 25% decrease 
of the final metallicity ZMWenv • 

It should be noted that the trend at high redshift is 
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Figure 9. MDF of the central galaxy calculated by GAMESH (solid 
black line) and in the IReion case (dotted red line) in the radial 
rang e 7 < r < 20 kpc. Black p oints indicate observed data taken 
from iBeers fe Christliebll2005l . 

certainly un-physical because the metallicity of the medium 
surrounding the galaxies on our box scale is expected to in¬ 
crease more progressively in redshift from very low Zmwbuv 
as correctly predicted by the IReion model. We thoroughly 
investigated what causes this trend, finding a concurrence of 
reasons. First, our N-body simulation does not provide a suf¬ 
ficient mass resolution to accurately predict the distribution 
of low mass mini-halos, which are expected to have played 
a relevant role in the early metal enrich ment of the MW 
and i t s dwarf satellites ( see for example Salvadori fc Ferraral 
120091 : ISales et al .11 20141 ; ISalvadori et al.l2014l ). As a conse¬ 
quence GAMESH tends to over-estimates the production/dis¬ 
persion of metals in the first star-forming objects. Sec¬ 
ond, as indicated by pure semi-analytic calculations that 
simultaneously reproduce d ifferent data-set s at z = 0 (see 
ISalvadori fe Ferraral (I 2 OI 2 I ); ISalvadori et al.l (I2014I )). a bet¬ 
ter description of both the stellar life-times and the in¬ 
homogeneous dispersion of metals into the MWenv is re¬ 
quired when the effect of an in-homogeneous reionisation 
are considered. Finally, on the radiative feedback side, the 
extension of the frequency range to a Lyman-Werner fre¬ 
quency band is necessary to correctly reproduce the star 
formation in LE-cooling halos. 

The Metallicity Distribution Function (MDF) of an¬ 
cient stars in the final MW-like galaxy, clearly reflects the 
different evolution of both the total SFR and ZMWenv at 
z > 8 — 10, where the majority of [Fe/H] < -3 stars are 
formed llde Bennassuti et al.ll2014 ). In Figure 9 we compare 
the MDFs of z = 0 stars in the radial range lkpc < r < 20 
kpc for the IReion and GAMESH runs. Black points with er- 
rorbars show the available data for Galactic halo stars in the 
same Galactocentric ranges (see also SF10). 

Both curves show a satisfactory agreement in the over¬ 


all trend of the observed data, and the differences be¬ 
tween the two runs are much smaller than the errors 
induced when considering different MW merger histories 
(ISalvadori. Schneider fe Ferraral 120071 ; Ide Bennassuti et al.l 
l2014h. While the GAMES H run fits the data for [Fe/H]> —3 
(see iBeers fe Christliebl 120051 ). it clearly overproduces the 
number of low-metallicity stars, [Fe/H]< —3. This is due to 
the slower progress of metal enrichment, as reflected by the 
shallow increase of ZMWenv in the redshift range 8 < z < 13 
(see the plateau in Figure 8). We interpret this issue as a 
clear indication that the current resolution of the simula¬ 
tion does not allow to accurately trace star formation and 
radiative feedback at the highest redshifts. 


5 CONCLUSIONS 

In this work we investigated the early formation process of a 
MW-like galaxy by using GAMESH, a novel tool combining for 
the first time in the literature, a N-body simulation with a 
detailed chemical and radiative feedback in a self-consistent 
framework. Many aspects of the formation process have been 
discussed along the cosmic time as the evolution of the SFR, 
the Pop III to Pop II transition, the statistics of progenitor 
halos in which star formation is suppressed, and finally the 
interplay between chemical and radiative feedback in setting 
up the chemical properties of the Local Group and the final 
MDF of the Galactic halo. 

Even with a low resolution N-body simulation, which 
lacks in statistics and does not allow to draw definite con¬ 
clusions on many aspects of the problem, GAMESH has been 
proven to be an ideal tool to provide both a comprehensive 
vision of the early formation of the Galaxy and its reion- 
isation process, and to include an accurate self-consistent 
treatment of chemical, mechanical and radiative feedback in 
numerical simulations. 

The capability of GAMESH to switch between pure semi- 
analytic treatment and inclusion of accurate RT enabled us 
to carefully compare and contrast all the results obtained 
by both models and to compare with observed properties 
of the Galaxy and the surrounding environment at z = 0. 
In addition, the GAMESH algorithm is very general and does 
not limit its field of applicability to a Milky Way formation 
process; it can be extended and adapted to a wide range of 
galaxy and star formation related topics. 

A future work, based on high mass and spatial resolu¬ 
tion simulation, is certainly needed to establish with accu¬ 
racy the epoch and extension of reionisation and its interplay 
with chemical feedback, the role of mini-halos and finally the 
signatures of global feedback processes on the stellar halo 
MDF and luminosity function of MW satellites. 
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